Primordial function and ambiguity in its 

determination 



Vladimir K. Petrov* 

N. N. Bogolyubov Institute for Theoretical Physics 
National Academy of Science of Ukraine 
252143 Kiev, Ukraine. 01.12.2003 



Abstract 

We study the possibility to reconstruct the primordial function for 
some periodic function. The procedure includes an analytical continuation 
of a discrete function for Fourier coefficients computation, that introduces 
an ambiguity. To establish conditions under which no ambiguity appears, 
the modification of Carlson theorem is suggested. In a case when no 
subsidiary condition is imposed, ambiguity in the definition of primordial 
function does not go beyond the functions which belong to space Q,' . 



1 Introduction 

In a previous work 1 we have put in correspondence to any tempered distribu- 
tion F {ip) some distribution 

00 

tF{^)^ ^ F{^ + 2T:n)^F{^) (1) 

n— — 00 

Recall that tempered distributions F {(p) G S' are defined as functional 

{F,<t>) = / F{^)<l){^)dp, (2) 



on a space S — {0} of fast decreasing test functions (f). Any F {(p) E S' can be 
represented as finite (m < 00) derivative F {tp) — -^p^f (</?) of some continuous 
tempered function ( \f {ip)\ < \ip\"' ; \ip\ ^ 00; a — const) Moreover, 
Fourier transform T [F {(p) , t] of any F {(p) G S' also belongs to space S" ^ |S] 
and every periodic distribution is tempered 7[. 

It is clear that S is linear operator, i.e. if E/ (ip) — f {(p) and Eg {(p) = g {p>), 
then £ [af (y) + (3g (^)] = af{^ ) + pgj^). 
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Procedure which we call the cyclization, is used not infrequently (see 
e.g. 0), but as a rule cychzation is applied to functions F (ip) that decrease 
quite rapidly with (p ±00 and F (ip) appears periodic. In ^] we considered 
the cases when this is not true and F {ip) ^ F {ip + 27r) . In this paper we handle 
only periodic F {ip) and consider the possibility to reconstruct function F (ip) 
(called here 'primordial') from the given F{ip), in other words, we consider 
procedure E~^F {(p) — F {(p) inverse to 1^. Since there is a family of functions 
uj (if) such that {ip) = 0, the inverse procedure is ambiguous and we estimate 
the degree of possible ambiguity. 

2 Primordial functions 

We will perform the cyclization procedure for F (p) represented as the Fourier 
integral hence an interchange of integration and summation order will be needed. 
Such interchange is justified in a case of uniform convergence. It is clear that if 
Ft £ S' , then the Abel-Poisson rcgularization for Fourier integrals 

/oo nOQ 
Fte'^'^dt^ / Fte'''^-'\'\dt (3a) 
-00 J —00 

and Fourier series^ 

00 00 
F{p)= F^e''^'^^ Yl -P'ne''^"""'"' (4) 



n— — 00 n— — 00 



provide a uniform convergence in and H3a|) . Abel-Poisson rcgularization also 
provides a standard representation (see e.g. |Q) for tempered distribution 

F{p) ^ F+{p + is)-F^{p-is); F{p) =F+{p + ie)-F^{p-ie) (5a) 
with 

nitoo 1 ±00 



F± {p) = / Fte'^'^dt; F± (p) = Tt:Fo ± V F„e'^" (6) 

Taking into account that 

00 00 



one can easily check that the apphcation of procedure to regularized integral 
lf5a|l gives 

poo 

„ ^ <y —00 



^Parameter s is a positive parameter, which is regarded as arbitrary small but finite and, 
if the opposite is not specified, it finally tends to zero. 
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The definition of product S {t — n) Ft for arbitrary Ft goes beyond the distri- 
bution theory and will be partly considered elsewhere. In this paper we consider 
only those Ft for which this product exists in the distribution theory framework. 
Recall that Dirac's (5-function is defined as 

/•CxD 

((5(t-s),(/)(t)) = / 6{t~ s)(j>{t)dt:^(f>{s) (9) 



for any (t) continuous in the small neighborhood of point t — s. 

After |S] we call the distribution Ft continuous in the area A (i.e. Ft G 
C° (A)) if there exists some function ft continuous in A and 

{{Ft -ft),cl, (t)) = j {Ft- ft) 4> it) dt = 0. (10) 

for any (j){t) E S with support located in A. As a A we choose the set of small 
neighborhoods (n — e,n + e) — A„.e of all integer n, where positive parameter 
e is arbitrary small, but finite. As it is shown in 0, expression 1)1 0|l defines 
ft uniquely almost everywhere. Further on Ft will be interpreted as piecewise 
'assembled' distribution which coincides with original Ft for t ^ ^n,e and we 
put Ft = ft for t e ^n,e- It is clear that in within the distribution theory such 
'assembled' distribution does not differ from the original one. 

A family of distributions {Ft} = S'^ continuous in the area A forms a linear 
space S'^ C S". Indeed, S'^ includes zero element and aiFi^t + oi2F2.t & S'^ for 
any Fj^^t £ S'^ and arbitrary complex numbers ak- For any Ft G S'^ integration 
in ^ may be done and we get 

F{ip)^ J2 / e^^'S{t-n)ftdt^ J2 (11) 



n— - 



In this case F {(p) can be represented as Fourier series (@J with coefficients given 
by the standard expression 

Fn^^J F {ip) exp {-iipn} dip (12) 

The simplest way to compute primordial function F {ip) is the straighforward 
extension of function Fn in (|12|l to the arbitrary value of parameter n —> t, i.e. 

— If''- 1 /■°° _ 

Ft^— F {p) exp {-ipt} dp ^ — F {p) e-xp {-ipt] dp) (13) 

27r J 



where 



We wish to stress the difference between Ft and Ft. Indeed, 

1 /'OO 1 p (2n4-l)7r 

Ft = — e-^^'F{p)dp = TT e-^^'F{p)dp 



m— — OO 

1 



27r 



e 



''^'F{p,t)dp (15) 
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where 

oo 

r(^,t)= e-^"''"F{v + 2n7T)^F{^) (16) 

771 — — oo 

This, in particular explains the difference in analytical properties of Ft and 
Ft. Expression shows that Ft is Fourier transform of distribution F (ip) 
with compact support, hence Ft is entire function. On the contrary, tempered 
distribution Ft may be singular. The only additional restriction is that Ft must 
be continuous on a small neighborhoods of all integer t — n. 

From (fl^ and ^ we see that EF {ip) = F (ip) so that F (ip) may be regarded 
as one of the primordial function for F{ip). Another primordial function for 
F {(p) may be obtained, if instead (|13|l one takes 

^'^'^i/T^^^^'""*'^^ ^^^^ 

where a is an arbitrary complex number. For all integer t = n we get ^ = Fn, 
but for arbitrary t the difference of extensions 

u;t = ~Ft = (e^'* - e"^^*) pt (18) 
does not turn into zero, if 

Pt^^ r F{ip-7T)e-''^'d<p^0 (19) 

Between F (ip) and F (ip) there exists a biunique correspondence. Despite 
the fact that no Fourier series reproduces F {ip) for \p\ > tt, expression (fT^ 
does not allow to discriminate coefficients F„ computed with F (cp) from the 
ones computed with truncated function F{p). However, extension F„ — * Ft 
with does not obligatory coincide with Ft given by p5|l . Indeed, from 
one can easily get 

— 1 f°° — ^ ^ sinTrfn — t) 

The difference between Fourier transforms of genuine and truncated functions.F 



cot = ^ r {F{p)-Fi^))e-''-'dp^Ft-Ft = Ft- £ f /"'"^" 

(21) 

generally turns into zero only for integer t. So the main problem in recon- 
struction of a primordial function is the estimation of admissible ambiguity in 
analytical continuation of 
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3 Uniqueness of straighforward extension 



If the class of functions is specified to which the required Ft belongs , this will 
impose restrictions on the analytical extension of Fn and may reduce ambiguity 
in the reconstruction of Ft from its values at integer t = n. It is known, that Fn 
may be uniquely extended on the whole complex plane , if Ft obeys the Carlson 
theorem^ conditions. According this theorem IHI, if i^t is an entire function and 
the indicator (called also the growth-measurihng function) defined as 

h {a; Ft) = lim^^oo^ In |Fc+roxp(m) | (22) 

obeys the conditions 

h{a;Ft)<X (23) 

for some A = const, and 

h[-^;Ft) +h{^; Ft) < 27:, (24) 

then 

F„=0; n = 0,l,2... (25) 

leads to Ft = 0. In other words, any entire function Ft is defined uniquely by 
its values at integer i = n = 0, 1, 2..., if Ft obeys conditions (|23|l and 124|) . 

If one takes the expression (|13|l for analytical extension Fn Ft, it will 
impose strong restrictions on Ft. Indeed, (jl3|l may be interpreted as Fourier 
transform of distribution with a support located in — tt < tp < tt, therefore 
Ft is an entire function and for \t\ oo 

h{a;Ft) <Tr\sina\ (26) 

so that condition (|23|l is satisfied. However, at the same time (|26|l does not 
exclude 

h[-^;Fn)+h[^:Fn)^27T, (27) 

hence condition (|24|) may be violated and Carlson theorem may fail to work. 
Incidentally, entire function 

1 f V— ^ r 1 \ r 1 J sin7r(t — n) 

ct = — / I anexp{iipn} j exp{-iipt\d(p =^ 2^ ^" n (t - n) 

\n— — oo / n— — OO ^ 

(28) 

obey H27(l and turns into zero for all t = 0,1, 2..., but it does not lead to Cf = 0. 

One can easily check that if F (if) turns into zero at some finite interval 
\ip — (po\ < e; < e < 27r, the analytical extension 113|l acquires the form 

= r ^F{^)e-^^'d^ (29) 



■^There are Carlson theorem alternatives (see e.g. |S] |5]) that differ by a set of such 
conditions. 
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and condition (|24|) is satisfied. Such case, however, is rather specific then rou- 
tine. 

Nonetheless, Carlson theorem may be easily modified for Fourier transforms 
of distributions with support located on interval of the length 2n. To be more 
specific, we consider F {ip) defined in (|14(l . Expression defines Fn uniquely. 
But if Fn is specified, series (@J uniquely defines F {ip) |H] and consequently 
F (if). Thereby if for all integer n (positive and negative) 

F„ = 0; 71 0,±1,±2.., (30) 

then from Q we get F {ip) = F {ip) — 0, which with (fT^ leads to Ft = 0. 



4 Linear space Q 

Here we estimate the maximum allowable ambiguity in computation of pri- 
mordial function when no subsidiary conditions are imposed. Any admissible 
extension Ft was to define Fourier coefficients Fn uniquely. This, in particular 
means that 

Fn = hm Ft = hm Ft (31) 

i.e. Ft was to be continuous in the area A defined above, so Ft is to belong 
to space S'^. Evidently the same will be true for the difference of any two 
extensions Fi^t — ^2,« = ^t- 

A family of distributions {w*} C S'^, such that Un — for all integer n forms 
linear space il. Indeed, includes zero element and aiuju +a2i-02,t G ^ for any 
£ ^ and arbitrary complex numbers at- 

It is clear that space f2 may include only those entire functions of exponential 
type which violate condition H24() . Moreover, the analytical extension of any 
LOn £ ^ with H13|) may give only trivial result ujt = 0, i.e. no cjj G il but = 
can be represented as (fT^ . 

Fourier transform 

/oc 
exp {i(pt} LOtdt (32) 
-oo 

relates space Q, = {w*} to linear space fi' — {lo {ip)}. Note, that if lo {(p) £ Vt' 
and tpQ is a constant which does not spoil the convergence of integral in (|32|l 
(e.g. Im(po = 0), then lo {tp + tpo) g . 
From {71) we obtain 

^ ^ />oo 

YiLo {ip) = uj{(p + 2'Kk)~ / exp {i(/3i} (5 (t — to) tJidt (33) 

Since ojt is continuous at integer t — n, integration in (|33f) may be done and, 
considering ujn = 0, we get 

oo 

f.uj {p) = t^„e''^" = w (v5) = (34) 
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i.e. procedure Q turns any of lo (ip) E fl' into zero. In other words, space 
fl' belongs to the kernel of the operator S. Note, that for any distribution 
t/i {(fi) such that JF~^ [tp {(p) , t] G S"^, the inverse proposition is true as well, i.e. 
(ip) = leads to tJj (ip) e Q'. 

It should be noted that procedure ^ turns any functions periodic with the 
period 2ttt into zero, if 1/t is noninteger number. Indeed, arbitrary function 
with the period r may be represented as 

oo 

$((^,t)= c„ (1/t) exp {mv^/r} (35) 

n— — c« 

SO, taking into account (Q we get 

oo oo 

E exp {icp/r} = exp {i {ip — 2iin) /r} — exp {iip/r^ 8 {m — 1/r) = 

n— — oo m— — oo 

(36) 

Now we were to summarize over all admissible r, i.e. over all noninteger t = 1/t. 
It means that c„ (l/r) = c„ (t) G f2 and consequently 

^ / c„ (t) exp (it = ^ C„ (.^) e 1)' (37) 



n— — oo n— — cxD 



and we see that no additional ambiguity is introduced. 

If l/r is integer, then Eexpjit/?} = exp {it/?} cothe/2 — > oo. Therefore, if 
some function / (</?) contains a term periodic with the period 27m, then E convert 
/ [p) into an infinite constant. Since we reconstruct primordial function F (p) 
from periodic function F (p) which does not contain infinite constants, such 
terms are out of consideration. 

For any lo (ip) e 57' one may define 

oo 

p{v)^p{^)+Pa{^); pa(v') =^w('^-^(2fc + l)) (38) 

fc=0 

where p{p) is arbitrary periodic distribution. With this definition one may 
write 

uj{p) = p{p + Tr) - p{p - tt) (39) 

that leads to^ 

LOt = (e^"* - e-"*) pt (40) 

We claim that singularity ol pt= T [p {p) , t] G S" at integer t = nis weaker 
than pole 

iOn = 2j lim (pf sinTri) = (41) 
that guarantees Wf G 51 and no additional restrictions on pt are imposed. 

•^In the particular case, considered above, uJt is given by 1181 and pt may be given by CHJ. 
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Note, that application of the procedure (^l to p ((p ± tt) with arbitrary p ((/p) 
does not obhgatory convert these functions into periodic ones 0. Therefore 
combination p(ip — t:) — p{ip + -k) does not obhgatory turn into zero after ap- 
phcation of the procedure iQ), but if condition H41I) is satisfied, it does. 

The change of e*'^* — e"*'^* in H4U|) for arbitrary polynomial in e*'^* and e~*'^* 
which turns into zero for all integer t leads to nothing but redefinition of pt- 

It is clear also that e.g. ujm.t = (e*'^* — e~*'^*) " gt G for any finite to, if 

Pt = (e"* - e-*'^*) qt obeys 

Let tempered distribution qt be bounded for the integers t — n and qt ^ 
at least for some integer t, i.e. qt ^ ft. Since 

hm Uit) , (e"* - e-"*)" = (42) 

for any cj) E S, then, by definition sequence f e il. However o^m.t converges 
to (7t ^ with TO — > oo. In other words, f2 is an incomplete space. 
Let us define the 'window' operator 51 as 

"'■'-{k lit: ■ 

which cuts out small neighborhoods A„^e of all integer t = n. It is natural to 
define 

/oo 
D,Ft exp {iipt] dt (44) 
-C30 

so for any F {ip) e 5' we get 



rn— — oo 



On the other hand, since 



with defined in (|10|) . one may write 



F{ip)^ / n - f]j Ftexp{i(/?t}dt = ^ / Ft exp {i(pt} dt 

(47) 

The continuity of Ft on A„_£ allows to apply the mean value theorem and rewrite 
ijTfll as 



1-?^)F(^) = 2^ f F„e-™.^^F(^). (48) 



7n— — oo 
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It is clear that the distribution [1 — 0.) F (w) vanishes with e ^ 0. Nonetheless, 



1-0) £^(..2.„). ^-y^^j;-'- . («, 



n— — oo 



hence, taking into account 

oo r ■ 1 

n— — OO 

(50) 

we find 

2 sin [{lp - 27rn) e] 
2- - 2.n) = ' 

n— — OO ^ ' 

and, consequently 

oo 

(l - 9}j SF ((p) = (^1 - f]) ^ F{^ + 2Tm) = F (^) (52) 

n— — oo 

Since 51^ = O, then both 57 and 1 — 57 may be considered as projectors. 
A family of tempered distribution S = with e S' , forms a linear 

subspace of ft' . For distributions Xt which are regular and continuous on all A„_e 
one can define a family of functions A = | ^1 — $7^ At|, which forms a linear 

space as well. Cyclization procedure ^ converts any A {(f) = T (1 — f2)At,(/J 



into a periodic function SA ((/?) = A ((/?) = A ((/? + 27r). 

Some linear transforms in spaces f7 and 57' are considered in Appendix. 

Now we turn to the LOt specified in (|21|) and try to compute the corresponding 
function pt. Taking into account 



X! / exp{-ja(t-n)}^ = ^ / (5 - ti^ exp {-iat} ^ = 1 

n— — OO n— — OO 

(53) 

one may rewrite H21|l as 

sinTTi ^ Ft- Fr, „ 

ujt = > — -1 (54) 

TT ^-^ t — n 

n— — oo 

that leads to 

P^ = ^Y1 -^-^{-IT (55) 

n— — oo 
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5 Conclusions 



In this paper we study a possibility to compute primordial function F ((^) for 
some periodic function F [ip) represented by Fourier series Q and presumably 
being a result of application of procedure ^ to F [ip). 

The reconstruction of primordial function demands an analytical continua- 
tion of Fourier coefficients, known only for integer i = n, to the whole real axis 
t. Such extension will introduce an ambiguity, unless the class of functions to 
which F [ip) belongs is not essentially restricted. 

As it is known, the extended function Ft is unique only if it obeys the 
conditions of Carlson theorem. In many cases such condition may appear too 
restrictive. We suggested a modification of Carlson theorem conditions for the 
case of extension given by ()13|l . which, as it seems to us, is the simplest way to 
compute primordial function F (</?) . 

In a case when no subsidiary condition is imposed, primordial function F (ip) 
may be reconstructed from F {(p) only up to arbitrary uj ((^) G fJ'. 



6 Appendix. Some linear transforms in spaces 
Q and Q' 

The distribution theory defines product or convolution only for narrow class 
of distributions. As it is stated in [S] any tempered infinitely differentiable 
function Kt is a multiplier in space S so there exists a product of such Kt with 
any tempered distribution. For some subspace of S' restrictions on Kt may 
appear to be weaker. 

Let us assume that for some family of ujf^^ G 51 there exists a product with 
some distribution Kt- It is natural to assume that KtJf^ e n. This wiU be 
true, indeed, if Kt is not 'too singular' at integers t — n, i.e. Kt and uJf^^ obey 
the condition 

hm Xtwp' = hm KtJ^^ = (56) 

i— >n+0 £— >n— 



for all integer n 

Exist 
forms T 



K] 



Existence of product KtLu\^^ means that convolution of the Fourier trans- 
it] (<^) s.ndT[Kt,p]^ K{^) 

OO 

K {a ~ Lp) J^^a) da = KJ^^ip) (57) 



exists. 

The main convolution property allows to write 



Sw^-^'l {ip) = w'^l {ip + 2Tm) = / K{a-ip) {a + 27rn) da 

-vn— ^— OO ^ — 

(58) 



n— ~oo 
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in other words 



(59) 



that also follows from Kt(^\^^ G fl- 
it is clear that a family of distributions fix — ^^[^^ | forms a linear subspace 

of the space ft. If in addition fix includes all KtLo].^^ then K converts il^ — 
{wl-^l i^)} into itself, i.e. i^f^^ = fl'j^. 
It is evident that in particular case 



1 



1 



1 



the convolution 



hJ"^ {if) = H 



wl^l (a) 



1 

—pv 

TT 



da. 



a — If 



(60) 



(61) 



may be interpreted as the Hilbert transform of cjl^l (ip). Symbol pv in (|61|l 
stands for principal value of the integral. 

From the fundamental property of the Hilbert transform we get 



(62) 



So the operator H converts subspace fl'^f C of functions cjl'^l which admit 
the Hilbert transform into itself. 
With the equation 



E 



(p — 27111 ± ie 



-exp{-n.} = ±-T — 



exp {±i (ip ± ie)} 



(63) 



one can directly check that T,H p] — 0. 

Recall that if cjl^l (ip) is 'ordinary' function, then to belong to fl'^j it must 
be piecewise continuously differentiable and there exists some constant A > 
such that for Ixl — + cx) 



(64) 



There are various approaches to Hilbert transform definition for distribu- 
tions. A brief review is given in ^U]. According to [3] the distributions / (a) 
which admit Hilbert transform H [f (a) , ip] may be defined as the functional 



{H[f (a),(p],0(^)) = 



1 

—pv 
n 



(a) 



: Ot-P 

{f{a),H[(p{p,),a]) 



da I {p>) dip 



(65) 
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on the test functions </) {ip) S which are infinitely differcntiable and cj) {(p) < 
0{l/\(p\) for \(p\ oo. Taking into account (|62l) we get H [(j)] S $ for any 
</> e $, or iJ [$] = <f>. In particular 



H 



(66) 



so H [a;(^)] S fi^. So taking into account (|62|l we see that Hilbert transform, 
indeed converts il'^j into itself. 

It should be noted that if condition (|64|l is not satisfied, Hilbert transform 
may be computed with some regularization. However in this case an important 
property of the Hilbert transform 

^- ruuf \ n - r H [f {a) ,ip - ie] - H [f (a) ,ip + ie] 

disc [H [f (a) , } = hm = / {(p) (67) 

may be broken. Indeed, Hilbert transform for / (a) = exp {iacr} 

H Te*"'", p] = -pv / da = e*°"P-pv / da = ie*"'^ signum (a) 

TT J_^a~p TT a 

(68) 

is the entire function of ip so disc {H [exp {iaa} ,p]} = 0. 

For some functions that does not obey H64|l the fundamental equation (j62|l 
may be broken too. Indeed, Hilbert transform of / = 1 is [1] = (see e.g. 
111115.2(1)). 

The reason is that we need the uniform convergence of the integral in H61|l to 
change the order of integration in (|62|l or the order integration and proceeding 
to limit in (|67|l . It is easy to prove that in the area < $ with arbitrary large 
but finite $ integral l|nT|l converge uniformly, if the condition is satisfied. 
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